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Abstract. The evolution of a binary star system by vari- 
ous analytic approximations of mass transfer is discussed, 
with particular attention payed to the stability of these 
processes against runaway on the thermal and dynami- 
cal timescales of the mass donating star. Mass transfer in 
red giant - neutron star binary systems is used as a spe- 
cific example of such mass transfer, and is investigated. 
Hjellming and Webbink's (1987) results on the dynamic 
timescale response of a convective star with a core to mass 
loss are applied, with new results. 

It is found that mass transfer is usually stable, so 
long as the the wind's specific angular momentum does 
not exceed the angular momentum per reduced mass of 
the system. This holds for both dynamical and thermal 
timescales. Those systems which are not stable will usu- 
ally transfer mass on the thermal timescale. Included are 
graphs illustrating the variation of §rj^- x ' = Cl with mass 
ratio in the binary, for various parameters in the non- 
conservative mass transfer, as well as evolutionary paths 
of interacting red giant neutron star binaries. 



as a screen in front of the other, as in EM Car, AO Cas, 
an d mos t notably t he bin ary PSR 1259- 63 (Jo hnston, et 
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1. Introduction 

The dominant feature in the evolution of stars in tight 
binaries, and the one which distinguishes it from the evo- 
lution of single stars is the presence of various forms of 
mass transfer between the two stars. Mass transfer occurs 
in many differe nt typ es of systems, to widely varying ef- 
fects (cf.: Shore 1994 ): contamination of the envelope of a 
less evolved star with chemically processed elements, as in 
Barium stars; winds from one star, which may be visible 



al, 1992; Kochanek, 1993; Melatos, et al., 1995; Johnston, 



et al., 1995 ); catastrophic mass transfer, by common enve- 
lope phase, as in W UMa systems; or a slow, steady mass 
transfer by Roche lobe overflow. 

Mass transfer is particularly interesting if one considers 
the evolution of a system with at least one degenerate 
star. In these cases, mass transfer produces spectacular 
effects, resulting in part from the intense magnetic and 
gravitational fields of the compact objects - pulsed X- 
ray emission, nuclear burning, novae outbursts, and so on. 
Also, since the mass transfer rates can be high, and orbital 
period measurements accurate, one may see the dynamical 
effects of mass transfer on the binary orbit, as in Cygnus 
X-3 (van Kerkwijk, et.al, |1992j; van den Heuvel |l994|). 
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In the case of cataclysmic variables (CVs) and low 
mass X-ray binaries (LMXBs), one has a highly evolved, 
compact star (CVs and LMXBs contain white dwarfs and 
neutron stars, respectively) and a less evolved main se- 
quence or red giant star. Mass transfer usually proceeds 
by accretion onto the compact object, and is secularly sta- 
ble. The transfer may be accompanied by a stellar wind 
from the mass-losing star, or ejection of matter from the 
accretor, as in novae and galactic jet sources. 

There are various unanswered questions in the evolu- 
tion of LMXBs into low mass binary pulsars (LMBPs, in 
which a millisecond radio pulsar is in a binary with a low 
mass white dwarf companion) and of CVs. Among these 
are the problem of the disparate birthrates of the LMXBs 
and the LMBPs (Kulkarni and Narayan 1988| ), estima- 
tion of the strength of X-ray heatin g induced winds from 
the d onor star (London, e t. al. |l981 ; London and Fl anner y 
1982t Tavani and London |l993| Banit and Shaham |l99^) , 
effects of X-ray heating on the red giant struc ture (Podsi- 
adlowski 1991 ; Harp az and Rappaport 1991 ; Frank, et.al. 
1992t Eamevxy, et.al. |1993[), etc. 



2 



G. E. Soberman, E. S. Phinney & E. P. J. van den Heuvel 



The problem of mass transfer in binaries by Roche 
lobe overflow has received a good deal of attention in the 
literature over the past few decades, typically in investi- 
gations of one aspect or another of orbital evolution or 
stability. Ruderman, et. al. (1989) examined mass trans- 



fer by isotropic winds and accretion, in investigating the 
evolution of short period binary X-ray sources with ex- 
treme mass ratios, such as PSR 1957+20 and 4U1820-30. 



King and Kolb ( |1995[ ) developed models with accretion 
and (typically) isotropic re-emission of transferred matter, 
in the context of the period gap in cataclysmic variables. 

The aim of this paper is to present a unified treat- 
mant of binary orbital evolution and stability, with rele- 
vant equations in a general form. The limits of pure modes 
of mass transfer and extreme mass ratios are presented 
and examined, to gain a qualitative understanding of slow, 
non-disruptive mass transfer in its various forms. 

The paper is organised as follows. Kinematic equa- 
tions for orbital evolution, assuming mass transfer through 
some well-specified mode(s), are derived in Sect. (^). Con- 
ditions for and regions of stability are derived in Sect. (|J). 
The theory is applied to LMXBs in Sect. (Ji|), and conclu- 
sions drawn in the final section. Two appendices contain 
further details on the unified models and on comparisons 
among the various limits of the models . 

2. Mass Transfer and the Evolution of Orbital Pa- 
rameters 

In this section, expressions for the variation of orbital pa- 
rameters with loss of mass from one of the stars are de- 
rived. In what follows, the two stars will be referred to as 
ra 1 and m 2 , with the latter the mass losing star. 

A binary composed of two stars with radii of gyration 
much less than the semimajor axis a will have an angular 
momentum: 



(1) 



where the period P is related to a and the total mass 
M T = m 1 + m 2 through Kepler's law: 



GM T 



(2) 



Here, the reduced mass y, = m 1 m 2 /M T and the mass 
ratio q = m^jm^. Tidal forces circularise the orbits of 
semidetatched binaries on timescales of 10 4 y <C T nuc p the 
nuclear timescale on which the binaries evolve (Verbunt 
and Phinney 1995). Furthermore, one expects stable mass 
transfer by RLOF to circularise orbits. Consequently, an 



1 The tidal timescale is also much shorter than the timescale 
of angular momentum loss by gravitational wave radiation, so 
binaries in which mass transfer is driven by angular momentum 
loss are also tidally locked. 



eccentricity e = will be assumed for the remainder of 
this work. 

Equations ([lj) and (Q) combine to give expressions for 
the semimajor axis and orbital period, in terms of the 
masses and angular momentum: 



L 



P 
2^ 



a = ( ~ \ {GM T )~ X 

3 

(GM T y 2 



with logarithmic derivatives: 

dlna = 2d\n{L/n)-d\nM T , 
dlnP = Zd\n{L/n)-2d\nM T 



(3) 
(4) 



(5) 
(6) 



There are various modes (or, borrowing from the ter- 
minology of nuclear physics, channels) of mass transfer 
associated with the different paths taken by, and destina- 
tions of mass lost from m 2 . 

The details of mass transfer must be considered, in or- 
der to calculate the orbital evolution appropriate to each 
mode. Several of these modes are described here, roughly 
following Section 2.3 of the review article by van den 



Heuvel (1994) 



The mode which is most often considered is accretion, 
in which matter from m 2 is deposited onto m 1 . The ac- 
creted component (slow mode) conserves total mass M T 
and orbital angular momentum L. In the case considered 
here, of a Roche lobe filling donor star, matter is lost 
from the vicinity of the donor star, through the inner La- 
grangian point, to the vicinity of the accretor, about which 
it arrives with a high specific angular momentum. If the 
mass transfer rate is sufficiently high, an accretion disk 
will form (c.f.: Frank, King, and Raine, 1985 ). The disk 
forms due to viscosity of the preferred fluid and trans- 
ports angular momentum away from and mass towards 
the accretor. As angular momentum is transported out- 
wards, the disk is expands to larger and larger circumstel- 
lar radii, until significant tides develop between the disk 
and the mass donor. These tides transfer angular momen- 
tum from the disk back to the binary and inhibit further 



disk growth (Lin and Papaloizou 1979). The mass transfer 
process is conservative if all mass lost from the donor (m 2 ) 
is accreted in this way. 

The second mode considered is Jeans's mode, which 



van den Heuvel calls, after Huang (1963), the fast mode; 



the third is isotropic re-emission. Jeans's mode is a spher- 
ically symmetric outflow from the donor star in the form 
of a fast wind. The best known example of Jeans's mode: 
the orbital evolution assoc iated with Type II supernovae 
in binaries (Blaauw, 1961 ) differs markedly from what is 
examined here. In that case, mass loss is instantaneous. 
A dynamically adiabatic outflow (m/m^> P orb ) is consid- 
ered here. 
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An interesting variant of Jeans 's mode is isotropic re- 
emission. This is a flow in which matter is transported 
from the donor star to the vicinity of the accretor, where 
it is then ejected as a fast, isotropic wind. The distinction 
between Jeans's mode and isotropic re-emission is impor- 
tant in considerations of angular momentum loss from the 
binary. Mass lost via the the Jeans mode carries the spe- 
cific angular momentum of the mass loser; isotropically 
re-emitted matter has the specific angular momentum of 
the accretor. 

Wind hydrodynamics are ignored. 

The fourth case considered is the intermediate mode, 
or mass loss to a ring. No mechanism for mass loss is hy- 
pothesised for this mode. The idea is simply that the ejecta 
has the high specific angular momentum characteristic of 
a circumbinary ring. 

The differences among each of these modes is in the 
variation of the quantities q, M T , and most importantly 
with mass loss. For example, given a q = 0.25 binary, 
mass lost into a circumbinary ring (a r > a) removes angu- 
lar momentum at least 25 times faster than mass lost by 
isotropic re-emission. This difference will have a great ef- 
fect on the evolution of the orbit, as well as on the stability 
of the mass transfer. 

In the following sections, formulae for the orbital evo- 
lution of a binary under different modes of mass loss are 
presented and integrated to give expressions for the change 
in orbital parameters with mass transfer. The different 
modes are summarised in Table (^J). 

Due to the formality of the following section and the 
cumbersome nature of the formulae derived, one should 
try to develop some kind of intuitition, and understand- 
ing of the results. The two sets of results (i.e.: one for a 
combination of winds, the other for ring formation) are 
examined in the limits of extreme mass ratios, compared 
to one another, and their differences reconciled in the first 
Appendix. 



Name of Mode Type 



Hhi osa /L 



Jeans/Fast Isotropic wind a 

— Isotropic Re-emission — f3 

Intermediate Ring 5 

Slow Accretion e 



A 

(1+9? 
g 

7 





Table 1. Names, brief description of, and dimensionless spe- 
cific angular momentum of the modes of mass loss explored in 
this paper. Parameters a, /3, 7, 8, e, and A are defined in Eqs. 
(7), (8), ||§, (37), and Q), respectively. The factor A in the 
fast mode is typically unity. 



2.1. Unified Wind Model 

Consider first a model of mass loss which includes isotropic 
wind, isotropic re-emission, and accretion, in varying 
strengths. Conservative mass transfer, as well as pure 
forms of the two isotropic winds described above, can then 
be recovered as limiting cases by adjusting two parameters 
of this model: 



P = - 



dTTlwind 

dm 2 
dm 2 



(7) 
(8) 



Here, dm W i n d is a small mass of wind from star 2 (isotropic 
wind), and dmi SO - r a small mass of wind from star 1 
(isotropic re-emission). Both amounts are expresses in 
terms of the total amount of mass dm 2 lost from star 2. 
A sign convention is chosen such that dm sta r is positive 
if the star's mass increases, and dm,fi ow is positive if re- 
moving matter (from m 2 ). Given the above, one can write 
formulae for the variation of all the masses in the problem: 



e m = l-a-0 
dlnm^ = —qe m dliim 2 

d\nM T = — — (a + f3)dlnm 2 
1 + q 

dlnq = (1 + e w q) d In m 2 



dlxifi = dlnm 2 



q 



1 



-d In q 



(9) 
(10) 

(11) 
(12) 

(13) 



Neglecting accretion onto the stars from the ISM and 
mass currents originating directly from the accretor, all 
transferred mass comes from the mass donating star (m 2 ) 
and a, (3, and the accreted fraction e ra all lie between and 
1, with the condition of mass conservation a + [3 + e m = 1, 
imposed by Eq. (9). For the remainder of this section, the 
subscript on e„ will be eschewed. 

If mass is lost isotropically from a nonrotating star, it 
carries no angular momentum in that star's rest frame. 
In the center of mass frame, orbital angular momentum 
L will be removed at a rate L = hrh, where m is the 
mass loss rate and h the specific angular momentum of 
the orbit. 

If the windy star is rotating, then it loses spin angular 
momentum S at a rate S = Ry^fl^rh, where is that 
star's rotation rate and i?^ the average of the square of 
the perpendicular radius at which the wind decouples from 
the star. Strongly magnetised stars with ionised winds will 
have a wind decoupling radius similar to their AlfVen ra- 
dius, and spin angular momentum may be removed at a 
substantially enhanced rate. 

The above consideration becomes important in the 
case of a tidally locked star, such as a Roche lobe filling 
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red giant. In this case, if magnetic braking removes spin 
angular momentum from the star at some enhanced rate, 
that star will begin to spin asynchronously to the orbit 
and the companion will establish tides to enforce corota- 
tion. These tides are of the same form as those between 
the oceans and the Moon, which are forcing the length of 
the day to tend towards that of the month. 

Strong spin-orbit coupling changes the evolution of or- 
bital angular momentum in two ways. First, for a given 
orbital frequency, there is an extra store of angular mo- 
mentum, due to the inertia of the star. Thus, for a given 
torque, L evolves more slowly, by a factor L/{L + S). The 
second is from the enhanced rate of loss of total angu- 
lar momentum, due to the loss of spin angular momen- 
tum from the star. Extreme values of Rw / R* allow the 
timescale for loss of spin angular momentum to be much 
less than that for orbital angular momentum. Thus, the 
second effect will either compensate for or dominate over 
the first, and there will be an overall increase in the torque 
due to this wind. 

This enhancement will be treated formally, by taking 
the angular momentum loss due to the fast mode to occur 
at a rate A times greater than what would be obtained 
neglecting the effects of the finite sized companion. 

Keeping the above discussion in mind, the angular mo- 
mentum lost from the system, due to winds is as follows: 



m 1 

M T 

JL 

Mo 



ego 



eg 

q 



1 + go 
q_ \ I I 



1 + eg 



eg 



1 + ego 
l + q J \l + eq 



(21) 
(22) 
(23) 



Furthermore, if the enhancement factor A is also 
constanta, then the angular momentum is an integrable 
function of q, as well: 



L 



1 



go 



1 



eg 



l + q J \l + eq 
where the exponents^ are given by: 



(24) 



B w = 
r — 



Aa , 
Aa + (3 
1 -e 
Aae 



1-e e(l-e) 



(25) 
(26) 

(27) 



Finally, substitution of Eqs. (22), (23), and (24) into 



Eqs. (3) and (4) gives expressions for the evolution of the 
semimajor axis and orbital period in terms of the changing 
mass ratio q: 



dL = {Aahi + [3hi)dm 2 , 



(14) 



which can be simplified by substituting in expressions for 
the hi\ 



Li/rrii 



M T — mi n 



hi 
hi 
dL 
<91nL 



q 2 L 
l + q m 2 
1 L 

l + q m 2 
Aa + (3q 2 L 
l + q m 2 
Aa + /3q 2 



dm 2 



i + q 



dim 



(15) 
(16) 
(17) 
(18) 
(19) 



So far, the equations have been completely general; a, 
(3, and A may be any functions of the orbital elements 
and stellar properties. Restricting these functions to cer- 
tain forms leads to simple integrable models of orbital evo- 
lution. If constant fractions of the transferred mass pass 
through each channel (a, (3 constant), then the masses are 
expressable as simple functions of the mass ratio q: 



in. 



q_ 



1 + ego 
1 + eg 



(20) 



«o 



P 



q_ 

go 
l + eq 



2A m -2 



1+gO 



1-2Z3„ 



1 



ego 
3 A, 



3+2C„ 



q_ 

go 
1 + eg 
1 + ego 



1 + go 



1-3Z3„ 



5+3C„ 



(28) 



(29) 



The derivatives of these functions may be evaluated 
either by logarithmic differentiation of the above expres- 
sions (Eqs. (28) and (29)), or by substitution of Eqs. (11), 
(12), (13), and (19) into Eqs. (5) and (6). Either way, the 
results are: 



d In a 
dlnq 



= 2{A W - I) + {I - 2B W )- 



-(3 + 2C U 



q 



l + eq' 



(30) 



2 This may not be the best approximation, if the orbit widens 
or narrows significantly and the torque is dominated by mag- 
netic braking. 

3 There is no problem when the denominator 1 — e — > 0, as 
the numerator vanishes at the same rate. For C w , the power 
laws become exponentials in the absence of accretion (e = 0). 
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<91nP 
dlnq 



3(A 



+(5 + 3C w ) 



9 



3B U 



1 + 9 



1 



"7 



(31) 



The reader should be convinced that the above equa- 
tions are correct. First, and by fiat, they combine to give 
Kepler's law. Second, they reduce to the familiar conser- 
vative results: 



a 
P 



1 



9 



i + 9o 

1 + 9 



1 



9o 



(32) 
(33) 



in the e — 1 limit. Finally, and again by construction, 
the formulae are composable. Ratios such as P/Pq are all 
of the functional form f(q)/f(qo), so if one forms e.g.: 
(P2/P0) = {P%/Pi){Px/Pq), the result is immediately in- 
dependent of the arbitrarily chosen intermediate point 

Pi = P(qi)- 

Note that the results for isotropic re- emiss ion obtained 
by Bhattacharya and van de n He uvel ( |1991 , Eq. (A. 6)), 
and again by van den Heuvel ( 1994 , Eq. (40)) are not com- 
posable, in the sense described above. The correct expres- 
sion for (a/ao) in the case of isotropic re-emission may be 
obtained by setting a = 0: 



Formulae describing orbital evolution can be obtained 
acording to the prescription of the previous section. If a 
fraction S of the mass lost from m 2 is used in the formation 
of a ring, then Eqs. (9) and (19) should be replaced as 
follows: 



e r = 1 - 6 
dL = 8h r dm 2 
— jSL/^idm 2 
= 7(5(1 + q)L/m 2 dm 2 
d\nL — 7<5(1 + q) dlnm 2 . 



(36) 



(37) 
(38) 



Following the procedure used in Sect. (2.1), and with 
5 = 1 (e r = 0): 



Mr 

L 

To 



a 
«o 

P 
Pb 



1 + 9 
l + 9o 



9o 



exp(7(q - q )) , 



q \ 2(7 - 1 Vl + 9 



</> ' / v 1 + qo 

3(7-1) 



1 + 9 
l + 9o 



exp(27(<j - q )) , 



exp(37(q- q Q )) 



(39) 
(40) 
(41) 
(42) 



a 



9o 
, 9 



l + qo\ / 1 + (!-/?)« 



1 



1 + (1 - 0)qo 



(34) 



Also no te tha t Eq. (A. 7) of Bhattacharya and van den 
Heuvel (1991) has its exponential in the denominator; it 
should be in the numerator. In van den Heuvel ( |1994|) , 
Eq. (37) should read d In J = ^ 1+ ^'°^ )g . In Eq. (38), 
an equals sign should repla ce the plus. (These corrections 
were also found by Tauris ( 1996 )). 

Thus, if there is mass loss from one star, with constant 
fractions of the mass going into isotropic winds from the 
donor and its companion, one can express the variation 
in the binary parameters, M T , a, and P, in terms of their 
initial values, the initial and final values of the ratio of 
masses, and these mass fractions. 

2.2. Formation of a Coplanar Ring 

Now consider a model in which mass is transferred by ac- 
cretion and ring formation, as described above. For con- 
creteness, follow a standard prescription (c.f.: van den 
Heuvel (1994)) and take the ring's radius, a r to be a con- 
stant multiple 7 2 of the binary semimajor axis. This ef- 
fectively sets the angular momentum of the ring material, 
since for a light ring, 



d In a 
<91ng 
dlnP 
dhiq 



2 7 (l + g)-2 + 
37(1 + 9) -3 + 



9 



1 + 9 
9 

1+9 



(43) 
(44) 



h r = L r /m r = y GM T a r = jL/fj, . 



(35) 



Unless the ring is sufficiently wide (a r /a ?Z 2), it will orbit 
in a rather uneven potential, with time-dependant tidal 
forces which are comparable to the central force. In such 
a potential, it would likely fragment, and could fall back 
upon the binary. Consequently, stability probably requires 
the ring to be at a radius of at least a few times a. A 'bare- 
minimum' for the ring radius is the radius of gyration of 
the outermost Lagrange point of the binary is between a 
and 1.25 a, dep ending on the mass ratio of the binary (Pen- 
nington, 1985 ). The ring should not sample the potential 
at this radius. For most of what follows, we will work with 
a slightly wider ring: a r /a — J 2 — 2.25 and 6 = 1. 

3. Linear Stability Analysis of the Mass Transfer 

Mass transfer will proceed on a timescale which depends 
critically on the changes in the radius of the donor star and 
that of its Roche lobe in response to the mass loss. The 
mass transfer might proceed on the timescale at which 
the mass transfer was initially driven (e.g.: nuclear, or 
orbital evolutionary), or at one of two much higher rates: 
dynamical and thermal. 
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If a star is perturbed by removal of mass, it will fall 
out of hydrostatic and thermal equilibria, which will be 
re-established on sound crossing (dynamical) and heat dif- 
fusion (Kelvin-Helmholtz, or thermal) timescales, respec- 
tively. As part of the process of returning to equilibrium, 
the star will either grow or shrink, first on the dynami- 
cal, and then on the (slower) thermal timescale. At the 
same time, the Roche lobe also grows or shrinks around 
the star in response to the mass loss. If after a transfer of 
a small amount of mass, the star's Roche lobe continues 
to enclose the star, then the mass transfer is stable, and 
proceeds on the original driving timescale. Otherwise, it is 
unstable and proceeds on the fastest unstable timescale. 

In stability analysis, one starts with the equilibrium 
situation and examines the small perturbations about it. 
In this case, the question is whether or not a star is con- 
tained by its Roche lobe. Thus, one studies the behaviour 
of the quantity 



AC 



m 5AR 
R 8m 



(45) 



which is the (dimcnsionlcss) variation in the difference in 
radius between the star and its Roche lobe, in response 
to change in that star's mass. Here AR is the difference 
between the stellar radius i?» and the volume-equivalent 
Roche radius r L . The star responds to this loss of mass on 
two widely separated different timescales, so this analysis 
must be performed on both of these timescales. 

The linear stability analysis then amounts to a com- 
parison of the exponents in power-law fits of radius to 
mass, R ~ m^: 



dlni?» 



d In to 
dlni?* 



d In to 
d In r L 



9 In to 



eq 



bin.evol. 



(46) 
(47) 
(48) 



where i?„ and to refer to the mass-losing, secondary star. 
Thus, i?» = r 2 and m = m 2 . Stability requires that after 
mass loss (8m 2 < 0) the star is still contained by its Roche 
lobe. Assuming Ai?2 = prior to mass loss, the stability 
condition then becomes 6AR 2 < 0, or £ t < (CCeJ- 11 
this is not satisfied, then mass transfer runs to the fastest, 
unstable timescale. 

Each of the exponents is evaluated in a manner consis- 
tent with the physical process involved. For £ s , chemical 
abundance and entropy profiles are assumed constant and 
mass is removed from the outside of the star. For £ e<n mass 
is still removed from the outside of the star, but the star 
is assumed to be in the thermal equilibrium state for the 
given chemical profile. In calculating Q L , derivatives are 
to be taken along the assumed evolutionary path of the 
binary system. 



In the following subsections these exponents are de- 
scribed a bit further and computed in the case of mass 
loss from a binary containing a neutron star and a Roche 
lobe-filling red giant. Such systems are thought to be the 
progenitors of the wide orbit, millisecond pulsar, helium 
white dwarf binaries. They are interesting, both by them- 
selves, and as a way of explaining the fossil data found 
in white dwarf - neutron star binaries. The problem has 
been treated by variou s aut hors, including Webbink, Rap- 
paport, and Savonije ( 1983 ), who evolved such systems in 
the case of conservative mass transfer from the red giant 
to the neutron star. 

3.1. Adiabatic Exponent: 

The adiabatic response of a star to mass loss has long been 
understood (see, for example, Webbink ( [1985 ) or Hjcllm- 
ing and Webbink (1987) for an overview), and on a sim- 
plistic level, is as follows. Stars with radiative envelopes 
(upper-main sequence stars) contract in response to mass 
loss, and stars with convective envelopes (lower-main se- 
quence and Hayashi track stars) expand in response to 
mass loss. The physics is as follows. 

A star with a radiative envelope has a positive entropy 
gradient near its surface. The density of the envelope ma- 
terial, if measured at a constant pressure, decreases as 
one samples the envelope at ever-increasing radii. Thus, 
upon loss of the outer portion of the envelope, the underly- 
ing material brought out of pressure equilibrium expands, 
without quite filling the region from which material was 
removed. The star contracts on its dynamical timescale, 
in response to mass loss. 

A star with a convective envelope has a nearly constant 
entropy profile, so the preceding analysis does not apply. 
Instead, the adiabatic response of a star with a convec- 
tive envelope is determined by the scalings among mass, 
radius, density, and pressure of the isentropic material. 
For most interesting cases, the star is both energetically 
bound, and expands in response to mass loss. 

Given the above physical arguements, the standard ex- 
planation of mass-transfer stability is as follows. A radia- 
tive star contracts with mass loss and a convective star ex- 
pands. If a convective star loses mass by Roche lobe over- 
flow, it will expand with possible instability if the Roche 
lobe does not expand fast enough. If a Roche lobe-filling 
radiative star loses mass, it will shrink inside its lobe (de- 
tach) and the mass transfer will be stable. 

This analysis is of only a meagre and unsatisfactory 
kind (Kelvin, 1894 ), as it treats stellar structure in only 
the most simplistic way: convective vs. radiative envelope. 

One can quantify the response of a convective star by 
adopting some analytic model for its structure, the sim- 
plest being an isentropic polytrope. This is a model in 
which the pressure P and density p vary as 



P(r) =Kp{r) 1+1 l , 



(49) 
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and the constituent gas has an adiabatic exponent related 
to the polytropic index through 7 = 1+1/n. Other slightly 
more realistic cases include those with 7^ 1 + 1/n, appli- 
cable to radiative stars; composite polytropes, with differ- 
ent polytropic indices for core and envelope; and centrally 
condensed polytropes, which are polytropes with a point 
mass at the center. These are al l considered in a paper 
by Hjellming and Webbink ( 1987 ). We use the condensed 
polytropes, as they are simple, fairly realistic models of 
red giant stars, which tend in the limit of low envelope 
mass to (pointlike) proto-white dwarfs which are the sec- 
ondaries in the low mass-binary pulsar systems. 

What follows is a brief treatmant of standard and con- 
densed polytropes, as applicable to the adiabatic response 
to mass loss. 

Scaling arguements give £ s for standard polytropes. 
Pressure is an energy density, and consequently scales as 
GM 2 R~ 4 . The density scales as MR~ 3 . The polytropic 
relation (Eq. j49[)) immediately gives the scaling between 
R and M. Since the material is isentropic, the variation 
of radius with mass loss is the same as that given by the 
radius-mass relation of stars along this sequence. Conse- 
quently, for polytropic stars of index n, 



71-1 



(50) 



In particular, £ 3 = —1/3 when n = 3/2 (7 = 5/3). □ 

The above approximation is fine towards the base of 
the red giant branch, where the helium core is only a small 
fraction of the star's mass. It becomes increasingly poor 
as the core makes up increasingly larger fractions of the 
star's mass, which happens when the star climbs the red 
giant branch or loses its envelope to RLOF. Mathemat- 
ically speaking, the scaling law that lead to the formula 
for ( s is broken by the presence of another dimensionless 
variable, the core mass fraction. 

A far better approximation to red giant structure, and 
only slightly more complex, is made by condensed poly- 
tropes, which model the helium core as a central point 
mass (see, e.g.: Hjellming and Webbink ( |1987[ )). Admit- 
tedly, this is a poor approximation, as concerns the core. 
However, the star's radius is much greater than that of 
the core, so this is a good first-order treatment. Further- 
more, differences between this approximation and the ac- 
tual structure occur primarily deep inside the star, while 
the star responds to mass loss primarily near the surface, 
where the fractional change in pressure is high. Overlap 
between the two effects is negligible. 

Analysis of the condensed polytropes requires integrat- 
ing the equation of stellar structure for isentropic matter 
(Lane-Emden equation), to get a function R(S, M c , M), 



4 This formula also reproduces the two following results. The 
gas giant planets (n = 1) all have approximately the same 
radius. Massive white dwarfs (n = 3) are dynamically unstable. 



and differentiating R at constant specific entropy S (adi- 
abatic requirement) and core mass M c (no nuclear evolu- 
tion over one sound crossing time). In general, the Lane- 
Emden equation is non-linear, and calculations must be 
performed numerically. The cases of n = and n = 1 are 
linear and quasi-analytic. The case of n = 1 is presented 
below, as a nontrivial, analytic example, both for under- 
standing, and because it can be used as a check of other 
numeric calculations of £„. 

The n = 1 Lane-Emden equation is (cf.: Clayton 



(1968)) 



dx dx 



(x) = -4>(x) . 



(51) 



Here, x — r y '2ttG / 'K = rjl is the scaled radial coor- 
dinate and <p = p/pc is the density, scaled to its central 
valuea. The substitution 4>{x) = x~ 1 u(x) allows solution 
by inspection: 



sxnjx) 

(x) = x e [0, 7rj 



(52) 



The polytropic radius is set by the position of the 
first root of <fi and is therefore at R = ir£. Similarly, 
M = 4ir 2 p c £ 3 . Eq. |[l]) shows that the density ((p) may 
be rescaled, without affecting the length scale, so R is in- 
dependant of M, and = Q. 

Alternately, the n = 1 polytrope admits a length scale, 
I which depends only on specific entropy, so the poly- 
trope's radius is independant of the mass. 

Generalising to condensed polytropes, Eq. ( p2[ ) sug- 
gests the extension: 



sin(:r + x ) _ rn 1 
<p{x) = x € [0, tt — x \ 



x 



(53) 



For this model, the stellar radius, stellar mass, and core 
mass are 



R 

M 



(tt - x )R , 
(tt - x )M , 
M(x -> 0) 
M sin(x ) , 



for some M and Rq. The core mass fraction 



M 



sin(z ) 

TT - Xq 



(54) 
(55) 

(56) 



(57) 



is a monotonic function in Xq, and increases from to 1 
as Xq increases from to tt. As the core mass fraction 
increases towards 1, the polytrope's radius decreases from 
ttRq to 0, so the more condensed stars are also the smaller 
ones. 



In general, p = pc(f) n . 
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The adiabatic R — M exponent C, a should be evaluated 
at constant core mass, as opposed to mass fraction. Thus, 
for condensed n — 1 polytropes, 



sin(xo) 



sin(xo) + (n — xq) cos(xq) 



(58) 



where xq is chosen to solve Eq. ( p7|) . This solution matches 
that given by Eq. ( |50| ) when there is no core. Furthermore, 
£ a is an increasing function of the core mass fraction, which 
diverges, as to = M c /M tends towards unity. These are 
general features of the condensed polytropes, and hold for 
polytropic indices n < 3. 

The procedure for calcu lating ( 3 is described in detail 
in Hjellming and Webbink (|1987| ). Results for a variety of 
core mass fractions of n = 3/2 polytropes are given, both 
m Table | and graphically, in Figs. ^ and |^. 



The function £ s (n 
by the function: 



3 . , 



3 1 



can be reasonably well fit 



(59) 



(Hjellming and Webbink, |1987| ), and to better than a per- 
cent by either of the functions (in order of increasing ac- 
curacy) : 



E 
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34.0 


0.11442 
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36.0 


0.09101 
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-0.18390 


38.0 


0.06925 


-1.0497 


-0.21759 


40.0 


0.04898 


-0.9785 


-0.24991 



Table 2. Adiabatic R — M relation vs. core mass fraction m c 



as in Hjellming and Webbink ( 1987 ), Table 3. Columns two 
and four are the core-mass fraction and mass-radius exponent, 
respectively. The parameter E in columns one and three is 
an alternate description of the degree of condensation of the 
polytrope, used by Hjellming and Webbink. The data presented 
here are in regions where the residuals of the fit formula Eq. 
(61) are in excess of 0.001. 



CsPHl 



1 / 1 — TO 



3 V 1 - mj 3V1 + 2to 



(sph = Csphi ~ 0.03m + 0.2 
as shown graphically in Fig. ||. 



1 + (1 - m)~ 



(60) 
(61) 



3.2. Thermal Equilibrium Exponent: 

In the case of a red giants, fits to data consistently show 
that the mass and luminosity depend almost entirely on 
the mass of the helium core (Refsdal and Weigert, |l970j ; 
cf.: Vcrbunt ( |l993|) ). Since the core mass changes on the 
nuclear timcscale, and we are interested in changes in the 
radius on the (much shorter) thermal timescale, the radius 
may be taken as fixed, giving = 0. 



3.3. Roche Radius Exponent: Q L 

Since the exponent Q L must be computed according to the 
evolution of the binary with mass transfer, it is sensitive 
to mass transfer mode, as are M Tl a, and P. The results 
for the various modes of nonconservative mass transfer 
are both interesting, and sometimes counterintuitive. It 
therefore makes sense to discuss them systematically and 
at some length. 

We rewrite d, in a form which depends explicitly on 
previously calculated quantities: 




0.4 0.6 
m = U/M 



Fig. 1. Plot of £ s versus core mass fraction, of an isentropic, 
n = 3/2 red giant star. As the fraction of mass in the core 
grows, the star becomes less like a standard polytrope. Im- 
portant to note is the crossing of the £ 3 and £ eq curves near 
m c /m s = 0.2. 
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-0.002 - 
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m - M /M g 



0.8 1 



d In r L 
d In m 2 

din a d\n(r L /a) d\nq 

d In m 2 d In q d In m 2 



(62) 



The derivatives and appear in Sect. (|), 

both for the unified model, and for the ring; as well as 
in a tabulated form in Appendix |A|. All that remains is 

9 * Vin L (/ a ' > ' wn i cn wm b e calculated using Eggleton's ( 1983 ) 
formula for the volume-equivalent Roche radius: 



r L /a 

din (r L /a) 
d\nq 



QA9q 



2/3 



0.6q 2 /3 + l n (i + g i/3) ' 

oV3 



1.2^/3 



1 



q i/3 o.6g 2 /3 + i n (i + g i/3) 



(63) 



(64) 



Fig. 2. 



Differences between various fit formulae and the 
The labeled curves are as fol- 
C), = 0.1(C S P H i - C«), and 



function Cs ( n — %'-> m 



lows: Zirw = 0.01 (Caw 




1 1.5 

q = m a /m, 



Fig. 3. Plots of ("s versus q, assuming a fiducial mi = 1.4Mq 
(q — m 2 /m 1 — m g /m x = (m c /m 1 )(m 2 /m c ) =(const)/(core 
mass fraction)). The three solid-line curves correspond, in as- 
cending order, to core masses of 0.14, 0.30, and O.45M0. One 
might choose to photoenlarge this figure, as well as (Q), et c, 
to use as overlays. 



To get an idea of how changing the mode of mass trans- 
fer effects stability, C, L has been plotted vs. q, for various 
models, in Figs. ||, [|, ||, and ^. 

One notices several things in these graphs. Figure ^ 
shows the extreme variation in C L (q) with changes in mode 
of mass transfer. Each of the three curves 'ring', 'wind', 
and 'iso-r' differs greatly from the conservative case. At 
least as important is the extent to which they differ from 
one another, based only on the way in which these three 
modes account for the variation of angular momentum 
with mass loss. When angular momentum is lost at an 
enhanced rate, as when it is lost to a ring or a wind from 
the less massive star (direct wind, at \ow-q; isotropic re- 
emission at high-g), the orbit quickly shrinks in response 
to mass loss and ( L is high. By contrast, in the case of the 
isotropic wind in the high-q limit, angular momentum is 
retained despite mass loss and the orbit stays wide, so ( L 
is lower than it is in the conservative case and the mass 
transfer is stabilised. 

Second, Fig. || shows that ring formation leads to 
rather high ( L , even for modest j 2 — a r /a — 2.25. The 
formation of a ring will usually lead to instability on the 
dynamical timescale ((„ < C, L ,ri ng )- The slower thermal 
timescale instability will occur only if the giant is rather 
evolved and has a high coremass fraction (with its high 

Delivery of mass to a ring would imply either a sub- 
stantially larger, simultaneous flow of mass through the 
stabler of the two Jeans's channels, or else orbital decay, 
leading to dynamical timescale mass transfer. 

Fig. ^ shows the variation of ( L due only to the dif- 
ferences between the isotropic wind and isotropic re- 
emission, in a family of curves with a + (3 — 0.8. At 
<7 = 1, stability is passed from the mostly isotropic re- 
emission modes at low q to mostly wind modes at high 



10 



G. E. Soberman, E. S. Phinney & E. P. J. van den Heuvel 



q. The q = 1 crossing of all these curves is an artifact of 
the model. Mass loss depends on the parameters in the 
combination a + /3. Angular momentum loss has a depen- 
dence Aa + (3q 2 . The two parameters a and (3 are then 
equivalent at q = V^4 = 1 . 

For <? < 1, mass transfer is stabilised by trading 
isotropic re-emission for wind, so families of a + /3 = const 
curves lie below their respective (/3 = 0, a) curves in this 
region of the q - Q L graph. This becomes interesting when 
one examines the /3 = curves, at high a. For total wind 
strengths of less than about 0.85 and q < 1, ( L < 0. Thus, 
for modest levels of accretion (at least ~ 15%), with the 
remainder of mass transfer in winds, red giant - neutron 
star mass transfer is stable on the thermal timescale, so 
long as q < 1. If the donor red giant has a modest mass 
core, so that £ s > 0, then the process will be stable on he 
dynamical timescale, as well. 

A family of curves with one wind of fixed strength and 
one wind of varying strength (as shown in Fig. ||, with (3 = 
0, and variable a), will also intersect at some value of q. It 
is easy to understand why this happens. For concreteness, 
take a constant-/? family of curves. Based on the linearity 
of the equations in a and (3, Cl(q~, ct, (3) = f(q; (3) + ag(q), 
with some / and g. If one evaluates Ct(?) &t a root of g(q), 
then the result is independent of a. 

Notice in each case, intersections of various unified 
wind model Cl{q) curves occur at q = O(l). One can 
understand this, in the framework of the arguments of 
Sect. (|A|). Depending on which side of q = 1 one lies, ei- 
ther m 1 or m 2 has the majority of the angular momentum. 
When q < 1, angular momentum will be lost predomi- 
nantly by the isotropic wind and £ L will be high. When 
q > 1, angular momentum losses will be modest and £ t 
low. The situation is reversed for isotropic re-emission. 
The two curves form a figure-X on the q - ( L diagram and 
must cross when niether q, nor l/q is too great; that is, 
near q = 1. 



: / 


/cons 


/ring 






/ / isor . 


Xwind / 






, , ,/ i / , , , 





0.5 



1.5 



2.5 



Fig. 4. £z, with all mass transfer through a single channel. 
For each curve, all mass from the donor star is transferred 
according to the indicated mode: conservative (cons) ; isotropic 
wind from donor star (wind); isotropic re-emission of matter, 
from vicinity of 'accreting' star (iso-r); and ring formation, with 
7 = 1.5. Since the unified model (winds+accretion) always has 
< a, (3, a + (3 < 1, the (, L (q) curves labeled cons, wind, and 
iso-r also form an envelope around all curves in the unified 
model. 



r = R Q expC^2c l y' 1 ) 
L = LoexpQ^^y 1 ) 



= Xe ( 



(65) 

(66) 
(67) 



4. An Example: Mass Transfer in Red Giant - Neu- 
tron Star Binaries 

Consider a binary system composed of a neutron star and 
a less evolved star. The less evolved star burns fuel, ex- 
panding and chemically evolving. If the binary orbital pe- 
riod is sufficiently short, the evolving star will eventually 
fill its Roche lobe and transfer mass to its companion. We 
now consider this problem, in the case where mass trans- 
fer starts while the donor is on the red giant branch (The 
so-called case B. See e.g.: Iben and Tutukov, 1985| . ) . 

The global properties of an isolated star are functions 
of the stellar mass and time from zero age main sequence. 
Alternately, red giant structure can be parameterised by 
total mass and core mass. Detailed models show that the 
dependence on total mass is weak, so the radius and lu- 
minosity are nearly functions of the core mass, only. We 



Here y = \xi{Am c /Mo), and the parameters a t , and c, 
are the result of fits to the red giant models (See Ta- 
ble g.). Equation (67) assumes that the red giant's lumi- 
nosity comes solely from shell hydrogen burning by the 
CNO cycle, which produces energy at a rate e CNO = 

10 18 c 



1983) 



5.987 • 10 i8 erg g" 1 (Webbink, et al. 

The red giant eventually fills its Roche lobe, trans- 
ferring mass to the other star. If mass transfer is stable, 
according to the criteria in Sect. |^, then one can manipu- 
late the time derivatives of r and r L to solve for the rate 
of mass loss from the giant. The stellar and Roche radii 
vary as: 



use the fit formulae by Webbink (1975), who writes: 



dr 
di 

dr L 



+ r C, — , 
m„ 



+ r L ( L 

m„ 



(68) 
(69) 
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Fig. 5. This illustrates the effect of an increasing fraction of 
the mass lost from the donor star (m 2 ) into a fast wind from 
that star. Solid lines correspond to a G 0.0(0.2)1.0 from top 
to bottom on the graph's right side; a — corresponds to 
conservative mass transfer. Note that at q ~ 0.72 , C, L {P = 0) 
is independent of q. 




2.5 



Fig. 6. The change in for all mass lost to isotropic 

winds, with varying fractions in isotropic wind and isotropi- 
cally re-emiiited wind. From top to bottom on the right side, 
the solid curves are for (a, f3) of (0.0,0.8), (0.2,0.6), (0.4,0.4), 
(0.6,0.2), and (0.8,0.0). Note that all curves cross at q = 1. 
Any set of curves a + f3 = const, will intersect at q = 1, as at 
this point, a and j3 have the same coefficients in ( L . 



Fig. 7. ( L for various a and /3, so that most transferred mass 
is ejected as winds. From top to bottom on the right side, 
the solid curves are for (a,P) of (0.6,0.2), (0.6,0.4), (0.8,0.0), 
(0.8,0.2), and (1.0,0.0). 



z 


ao 


ai 


ai 


a.3 


0.02 
10" 4 


3.50 
3.27 


8.11 
5.15 


-0.61 
4.03 


-2.13 
-7.06 


Z 


CO 


Cl 


C2 


C3 


0.02 

io- 4 


2.53 
2.02 


5.10 
2.94 


-0.05 
2.39 


-1.71 
-3.89 



Table 3. Parameters fitted to a series of red giant models, both 
for Pop I (Z = 0.02, X = 0.7) and Pop 11 {Z = 10~ 4 X = 0.7). 
Data are transcribed from Webbink, el. ai, (1983), and are 
applicable over the ranges y € (—0.4,0.6) and y £ (—0.2,0.4) 
for Po p I and Pop II, respectively. Pop I fi gures due to Webbink 
jl975|) ; Pop II from Sweigart and Gross ( |i978| ). 



So long as the star remains in contact with its Roche lobe, 
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m g 
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r L 
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r 


dt 





d In r dim 



dt 



dr L 
dt 



(70) 



(71) 



The first term in Eq. (69) takes into account changes in 
the Roche radius not due to mass transfer, such as tidal 
locking of a diffuse star or orbital decay by gravitational 
wave radiation. For the models considered here, the Roche 
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lobe evolves only due to mass transfer. Equations (69), 
(70), and (71) then reduce to the two equations: 



dint 
~dt 



The additional relation: 
d In r d In r d In m r 



dt 



d In m c dt 



(72) 
(73) 

(74) 



is a consequence of the star being a red giant. 

The equations for the evolution of the core mass, red 
giant mass, neutron star mass, orbital period, and semima- 
jor axis (obtainable from Eqs. (67), (73), @, (11), (31), 
and (30)), form a complete system of first order differ- 
ential equations, governing the evolution of the red giant 
and the binary. The core mass grows, as burned hydrogen 
is added from above. This causes the star's radius to in- 
crease, which forces mass transfer, increasing a and P, at 



of the semi-detatched requirement r 2 = r L = a ■ (r L /a). 
Plots of various quantities versus red giant core mass show 
the evolutionary histories of binary systems, in Figs, pi ^, 
and [l(]. Population I stars were used in all calculations 
(see Table %Z = 0.02). 

Another interesting case of mass transfer is isotropic 
re-emission at the minimum level necessary to ensure Ed- 
dington limited accretion: 



(3 = max(0, fn g /rh > 



1) 



(75) 



Re-emission would presumably be in the form of pro- 
peller ejecta (Ghosh and Lamb, ( 1978 ) or a bipolar out- 
flow, such as the jets seen in the galactic superluminal 
sources. The evoution of a binary, subject to this con- 
straint, is shown in Fig. [H]. Comparing Fig. [lljwith Figs. 
H and ^, one sees that the low core mass systems, with 
their corresponding low mass transfer rates, mimic conser- 
vative systems. The faster evolving systems behave more 
like systems with pure isotropic re-emission. 

Almost all ring-forming systems are unstable to ther- 
mal and/or dynamical timescale runaway of mass transfer, 
and are not displayed. 



4-1- The Code 



5. Conclusion 



The program used to generate the numbers presented here 
follows the prescription outlined by Webbink, et al. 1983| . 
in the treatment of initial values, numeric integration of 
the evolution, and prescription for termination of mass 
transfer. T he so lar mass, radius, and luminosity are taken 
from Stix (|1991|) . 

Initial values were provided for the masses of the neu- 
tron star, the red giant, and its core at the start of the 
contact phase (red giant filling its Roche lobe): m x , m g , 
and to c , respectively, as well as for the parameters (a, j3, 
7 ...) of the mass transfer model. A tidally- locked system 
was assumed and spin angular momentum of the stars ne- 
glected. The program then solved for the orbital period P 
and separation, a, using relations (||), (63), and (65). 

Integration of Eqs. (10), (6), (67), and (73) was 
performed numerically by a fourth-order Runge-Kutta 
scheme (cf.: Press, et al. 1986| ) with time steps limited by 
Am g < 0.003m g , Am c < 0.001m c , and Am e < 0.003m e . 
The first two criteria are those used by WRS; the last was 
added to this code, to insure that care is taken when the 
envelope mass, to c , becomes small toward the end of the 
integration. 



Detailed numeric calculations (Taam, 1983) show that 
a red giant cannot support its envelope if m c /m g < 0.02. 
The code described here follows that of WRS and termi- 
nates mass transfer at this point. 

An overdetermined system of m c , m gl m x , and P was 
integrated numerically by the code, allowing for consis- 
tancy checks of the program. Tests performed at the end of 
the evolution included tests of Kepler's law (Eq. (||)), an- 
gular momentum evolution (Eq. (29) or (42)), and a check 



Previous works in the field of mass transfer in close bina- 
ries usually centered on conservative mass transfer, even 
when observation directs us to consider mass loss from 
the system, as is the case with white dwarf neutron star 
binaries. If these systems form from mass transfer in red 
giant neutron star systems, then one must find a way to 
start with a secondary star sufficiently massive to evolve 
off the main sequence in a Hubble time, (m 2 > 1M®), in 
a binary with a m x > 1.3M Q neutron star and reduce the 
secondary's mass by ~ 0.6M Q , while keeping the neutron 
star below ~ 1.45M . To do this, some 3/4 of the mass 
from the secondary must be ejected from the system. Re- 
sults here indicate that it is possible to remove this much 
mass in winds, while maintaining stable mass transfer on 
the nuclear timescale. 

For the most part, nonconservative mass transfer, in 
which mass is lost in fast winds, mimics the conservative 
case. For the typical initial mass ratios (q « 0.7), ( L ranges 
from about -0.7 to -0.2. The rate of mass transfer is given 
by equation (73), and is longest at the start of mass trans- 
fer, when m c is low. The total time is therefore set almost 
entirely by m e (0) and Cl(0), and so differs by maybe a 
factor of 3, over all possible wind models. 

It is worth repeating that changes in M T arising from 
nonconservative evolution do not alter the relationship be- 
tween white dwarf mass and binary period by more than 
a couple of percent. This is true, by the following argue- 
ment. The final red giant and core masses differ only at 
the few percent level. Approximately, then, the red giant 
mass sets the red giant radius and therefore the Roche 
radius. In the (low q) approximation used by Paczyriski 
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(1971), orbital period is a function of r L and m 2 , only. P 
vs. m c is a function of the final state, alone. Using Eggle- 
ton's formula instead of Paczynski's introduces only a very 
weak dependence on the mass of the other (neutron) star. 
In the end, the theoretical motivation for the existance of 
a P - m c relation is significantly more solid than, say, our 
knowledge of the red giant R - m c relation, on which the 
exact P - m c curve depends. 

The exact mode of mass transfer will effect P/Po, as 
is evident from Eq. (29). This could be important, in sta- 
tistical studies of white dwarf neutron star binaries, and 
trying to predict the distribution of P from the initial mass 
function, and distribution of initial orbital periods. This is 
dependant, of course, on the development of a quantitative 
understanding of the common envelope phase. 

Finally, and probably the most useful thing, is that if 
one assumes only accretion and wind-like mass transfer, 
then most every binary in which mass is transferred from 
the less massive star is stable on both dynamical and ther- 
mal timescales. If the mass donor has a radiative envelope 
(not treated here), it will shrink in response to mass loss, 
and lose mass stably. If the donor has a convective enve- 
lope, a modestly sized core will stabilize it sufficiently to 
prevent mass loss on the dynamical timescale. Only if one 
has a very low mass (or no) core, will the mass transfer be 
unstable on the dynamical timescale, and then, only for 
a — ► 1. High values of a, and low core mass in a convec- 
tive star may lead to instability on the thermal timescale, 
if the mass ratio, q is sufficiently low (see Fig. |J). 



Acknowledgements. This work was supported by NSF grant 
AST93- 15455 and NASA grant NAG 5-2756. 

References 

Banit, M. and Shaham, J. 1992, ApJ, 388 L19 
Bhattacharya, D. and van den Heuvel, E.P.J. 1991, Phys. Rep. 
203, 1 

Blaauw, A. 1961, Bull. Astron. Inst. Neth., 15, 265 
Clayton, D.D. 1968, Principles of Stellar Evolution and Nucle- 
osynthesis, University of Chicago Press. 
Eggleton, P.P. 1983, ApJ, 268, 368 

Flannery, B.P. and van den Heuvel, E.P.J. 1975, A&A, 39, 61 
Frank, J., King, A.R., and Lasota, J. P. 1992, ApJ, 385, L45 
Frank, J., King, A.R., and Raine, D.J. 1985, Accretion Power 

in Astrophysics, Cambridge University Press. 
Ghosh, P. and Lamb, F. K. 1978, ApJ, 223, L83 
Hameury, J.M., King, A.R., Lasota, J. P., and Raison, F. 1993, 

A&A, 277, 81 
Harpaz, A. and Rappaport, S. 1991, ApJ, 383, 739 
Hjellming, M.S. and Webbink, RF. 1987, ApJ, 318, 794 
Huang, S.S. 1963, ApJ, 138, 342 

Iben, I. Jr. and Tutukov, A.V. 1985, ApJ, Supp. Ser., 58, 661 
Johnston, S., Manchester, R.N., Lyne, A.G., Bailes, M., Kaspi, 

V.M., Guojun, Q., and D'Amico, N. 1992, ApJ, 387, L37 
Johnston, S., Manchester, R.N., Lyne, A.G., DAmico, N., 

Bailes, Gaensler, B.M., andNicastro, L. 1995, MNRAS, 279, 

1026 



Kelvin, Sir William Thomson, Lord. 1891-1894, Popular Lec- 
tures, Macmillan & Co. 
King, A.R. and Kolb, U. 1995 ApJ, 439, 330 
Kochanek, C.S. 1993, ApJ 406, 638 
Kulkarni, S.R. and Narayan, R. 1988, ApJ 335, 755 
Landau, L.D. and Lifshitz, E.M. 1951, The Classical Theory 

of Fields, Addison- Wesley Press, Inc. 
Lin, D.N.C. and Papaloizou, J. 1979, MNRAS, 186, 799 
London, R., McCray, R., and Auer, L.H. 1981, ApJ, 243, 970 
London, RA. and Flannery, B.P. 1982, ApJ, 258, 260 
Melatos, A., Johnston, S., and Melrose, D.B. 1995, MNRAS 
275, 381 

Paczyriski, B. 1971, Ann. Rev. A&A, ARAA 9,183 
Pennington, R. 1985, in Interacting Binary Stars, Pringle, J. 

E. and R. A. Wade, eds., pl97. 
Podsiadlowski, Ph. 1991 Nature, 350, 136 

Press, W.H., Flannery, B.P., Teukolsky, S.A., and Vetter- 
ling, W.T. 1986, Numerical Recipes, Cambridge University 
Press. 

Refsdal, S. and Weigart, A. 1970 A&A, 6, 426 
Ruderman, M., Shaham, J., and Tavani, M. 1989, ApJ, 336, 
507 

Shore, S.N. 1994, in Interacting Binaries (Saas-Fee 22), Shore, 

S.N., et al., eds., pi 
Stix, M. 1991, The Sun, Springer Verlag. 
Sweigart, A.V. and Gross, P.G. 1978, ApJ Supp. Series, 36, 

405 

Taam, R.E. 1983, Ap. J, 270, 694 

Tauris, T., 1996, A&A, 315, 453. 

Tavani, M. 1991, ApJ 366, L27 

Tavani, M. and London, R. 1993, ApJ, 410, 281 

van den Heuvel, E.P.J. 1994, in Interacting Binaries (Saas-Fee 

22), Shore, S.N., et al., eds., p263 
van Kerkwijk, M.H., Charles, PA., Geballe, T.R., King, D.L., 

Miley, G.K., Molnar, L.A., and van den Heuvel, E.P.J. 1992, 

Nature, 355, 703 
Verbunt, F. 1993, Ann. Rev. A&A, 31, 93 
Verbunt, F. and Phinney, E.S. 1995, A&A, 296, 709 
Webbink, RF. 1975 MNRAS, 171, 555 

Webbink, RF. 1985, in Interacting Binary Stars, Pringle, J.E. 

and R.A. Wade, eds., p39. 
Webbink, R.F., Rappaport, S., and Savonije, G.J. 1983, ApJ, 

270, 678 



A. Mass Transfer Models With Extreme Mass Ra- 
tios 



The results of Sees. |2.1| and |2.2| are examined here in the 
limiting cases of q tending toward zero and infinity, to 
explore the relations between the different modes, and the 
connections to well-known toy models. This is done to 
develop an intuition which may be used in comparing the 
stability of mass transfer by the various modes. 

In considering extreme values of the mass ratio, one 
makes the reduced mass approximation, regarding all the 
mass as residing in one star of fixed position and all the 
angular momentum in the other, orbiting star. Errors are 
only of order q (or 1/q, if this is small). One might keep the 
Solar System in mind as a concrete example. The Sun has 
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all the mass, and the total angular momentum is well ap- 
proximated by Jupiter's orbital angular momentum. The 
total mass M T ~ M© and the reduced mass /i ~ Mj, with 
fractional errors ~ Mj/Mq ~ 10 -3 . 

Retention of the factors A (Eqs. (19), et seq. of 
Sect. O) and 5 (Eq. (38), of Sect. |2~^) allow for broader 



comparisons between the two models. Formulae appropri- 
ate to the two extreme limits are in Table |]. Columns 2 
and 3 pertain to the wind models; 4 and 5 to the ring. 

It is now straightforward to compare the two models 
of mass transfer and loss in the extreme limits where they 
should be equivalent. Comparisons will be made first in 
the test-mass is donor limit, where the unified and ring 
models correspond exactly. Second is a treatment of the 
opposite limit, where the two models give different but 
reconcilable results. 

First, the winds and ring models are equivalent in the 
q = limit, upon identification of a with S and A with 7H. 
This symmetry is not difficult to understand. In the q = 
limit, only the direct isotropic wind and ring remove an- 
gular momentum. In this case, isotropic re-emission is an 
isotropic wind from a stationary source, and accretion is 
always conservative of both mass and angular momentum. 
In each torquing case (direct wind and ring), the ejected 
mass removes specific angular momentum at an enhanced 
rate — A in the winds model, 7 in the ring model. Thus, 
it makes sense that wherever A and 7 appear in the equa- 
tions, they are in the products Aa and 7$, the rate of 
angular momentum loss per unit mass lost from m 2 . 

More interesting is that, in the strict q = limit, the 
parameters a and 5 are found only in the combinations 
Aa and jS. The independance from other combinations 
of parameters can be understood by examining the ratio 



V 



5 In M T 



9 In (L/fi) 



(Al) 



which tells the relative importance of angular momentum 
and mass losses in the orbit's evolution. For small donor 
mass, i] is small and mass loss without angular momen- 
tum loss is unimportant. The strict mass loss term will 
be important only when the coefficient of the d In (L/fx) 
term, Aa — 1 or j5 — 1, is of order q or less. This happens, 
for example, in the Jeans' mode of mass loss (q <C 1, 
A = a = 1), discussed below. 

Even neglecting the questions of stability important 
to tidally induced mass transfer, the situation is differ- 
ent when the mass-losing star has all the mass and al- 
most no angular momentum. In this case, both mass and 
angular momentum loss are important. Since the more 
massive body is the mass donor, a non-negligible fraction 
of the total mass of the system may be ejected. Further- 
more, the angular momentum per reduced mass changes 
via isotropic re-emission and ring formation. Therefore, 



both mass and angular momentum loss play a role in the 
dynamics. The test of relative importance is 77, which in 
the limit of large q, goes as 1/q. Again, in the case of ex- 
treme mass ratios, changes in L/fi dominate the course of 
evolution. As before, there are times when the coefficient 
of the i9 In (L/fi) term, 1 — a or jS + (1 — S), is of order 
q, or less. In this case, above arguments fail and the strict 
mass loss term ainM T must be included. 

The evolution of the angular momentum per reduced 
mass shows a significant difference the winds and ring 
models in their g > 1 limit. In the winds case, one may 
write q(l — a) = q(l — a — (3) + q/3. The first part is the 
fractional rate of increase of /i, with loss of mass from the 
donor star. The second is the rate at which angular mo- 
mentum is lost from the system, by isotropic re-emission. 
Since g> 1, the mass losing star is nearly stationary and 
mass lost in a direct wind removes no angular momentum. 

The ring + accretion case is slightly different, but 
may be interpreted similarly. The term q(l — S) replaces 
q{l — a — P) as the fraction of mass accreted onto the 
first star; i.e.: ain/i/ainm 2 . The term q(3 is replaced by 
q"fS, as the rate of loss of angular momentum. Although 
the donor star is stationary in this 9 > 1 ring case, an- 
gular momentum is still lost. This is due to the particular 
construction of the ring model, in which the angular mo- 
mentum removed is proportional to the system's Lj /i, not 
to the specific angular momentum of either body. Thus, 
in the ring model, even as hi tends towards zero with 
increasing q, ejected mass will carry away angular mo- 
mentum. This is the final caveat: the ring model is just a 
mechanism for the rapid removal of angular momentum 
from a binary system, and one should keep this in mind, 
particularly when q 1. 

One might also notice for both models, that if angular 
momentum loss is not overly efficient (^4 and 7 not too 
much greater than 1), then mass loss from the test mass 
widens the orbit and mass loss from the more massive star 
shrinks the orbit, just like in conservative mass transfer. 
Also, 77 is not everywhere small. In particular, when q is 
of order unity, 77 is also and the d In M T term becomes as 
important as d In (L/fi) in the equations of motion. 

A.l. Jeans's Modes 

Often, one talks of the Jeans's mode of mass loss from a 
binary, in which there is a fast, sperically symmetric loss 
of mass. The Jeans' mode has two limits. One is a catas- 
trophic and instantaneous loss of mass, as in a supernova 



Actually, one must only identify the products of mass frac- 
tion lost and relative specific angular momenta: Aa and 7<5. 



event (see van den Heuvel (1994), for a discussion). In this 
case, if the orbit is initially circular and more than half 
the total mass is lost in the explosion, the system unbinds. 
This is can be explained, based on energetics. Initially, the 
system is virialised with < T >= —1/2 < V >. The loss of 
mass does not change the orbital velocities, so the kinetic 
energy per unit mass remains the same. The potential en- 
ergy per unit mass is proportional to the total mass, so 
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if more than half the mass is lost, E = T + V > and 
the orbit is unbound. A more detailed analysis may be 
done, and will give ratios of initial to final orbital periods 
and semi-major axes, for given initial to final mass ratios 
(Blaauw ( |1961 ); Flannery and van den Heuvel ( 1975| )). 

The equations in this paper will not give the 'stan- 
dard' Jeans solution and unbound orbits. Unbinding the 
orbit from elliptical to hyperbolic requires e > 1, where 
e = has been assumed from the outset. The argument 
used in the preceding paragraph seemingly necessitates 
the unbinding of the orbit with sufficient mass loss, but it 
is not applicable here, as it also assumes a conservation of 
mechanical energy. Mechanical energy is not conserved in 
the above calculations (section (2T), for example), as the 
presence of dissipative forces to damp e to have been 
assumed. 

The other limit of Jeans's mode is mass loss by a fast 
wind, on a timescale slow compared to the orbital pe- 
riod. In this case, there is no preferred orientation for the 
Runge-Lenz vector (direction of semimajor axis in an ec- 
centric orbit), and the orbit remains circular throughout 
the mass loss, with M T a — const. This is the limit of 
the Jeans mode which our calculations reproduce. Jeans's 
mode is an example of a degenerate case (mentioned 
above; here, Act — 1 ), where <91n(L//i) = 0. In this 
case, one may simply apply equation (3), and see that 
GM T a = (L/ '/i) 2 = const. 

B. Extensions to models considered above 

In the interest of completeness, a five-parameter model 
of mass transfer, combining winds from both stars, ring 
formation, a nd a ccretion is presented. The treatment in 
Sects. H and 3J3 is followed. Modifications necessary for 
inclusion of other forms of angular momentum loss, such 
as Lg-Rj due to gravitational radiation reaction, are also 
discussed. 

Construction of this model is straightforward, and re- 
sults from the inclusion of the various sinks of mass and 
angular momentum, due to various processes. The noncon- 
servative part of each model makes its own contribution 
to the logarithmic derivatives of M T and L: 



e = 1 — a - 
d In M x = (1-e) 



<1 



1 



q 



s, 

- d In m 2 , 

+ 7<5(1 + q) ) <91nm 



(Bl) 
(B2) 

(B3) 



where each of the variables retains its old meaning. Re- 
placing the old definition of e in Eq. (9) with the new 
definition in Eq. (Bl) makes Eqs. (10) - (13) applicable to 
this model, as well. 

Contributions from other evolutionary processes, such as 
gravitational wave radiation reaction; realistically pre- 
scribed stellar winds simple prescription) ; tidal evolution, 



et c. can also be added. Each will make its own contribu- 
tion to the angular momentum and total mass loss. For 
example, orbital decay by gravitational radiation reaction 



(Landau & Lifshitz, 1951) can be included as an other sink 



of angular momentum: 



dt GR 



5c 5 



(GM t ) 6 (L/m)" 



(B4) 



In most cases, this type of physics can be modelled 
as an intrinsic r L (the second term in Eq. (69)). Like the 
intrinsic stellar expansion term of Eq. (68), this kind of 
evolution occurs even in the absence of mass transfer. 
Therefore, the convenient change of variables from t to 
q used in Sect. (Q) introduces singularities when evolution 
takes place in the absence of mass transfer. The equa- 
tions of Sect. (Q) still hold, but only for those phases of 
the binary's evolution durring which tidally-driven mass 
transfer takes place. 

We temporarily neglect these complications and con- 
sider mass transfer via isotropic wind, isotropic re- 
emission, and formation of a ring, with mass fractions a, 
/3, and 5 respectively. The remainder of the mass transfer 
(the fraction e = 1 — a — f3 — S) goes into accretion. The 
ratio h r /hbi n is 7, where j 2 — a r /a. Note that 5, 7, 0, 
and a are all used as before; e should still be regarded as 
the accreted fraction. 
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Taking a, (3 = 0, gives formulae for a ring of stren; 



(B7) 
(B8) 
(B9) 

(BIO) 

rth 6: 



«o 



P 



2(1-A r 



1 



1 



! + (!-% 
1 + (1 - 5)q a 
3(1-A r 



go, 

3+2C r 



(Bll) 



i + g 



16 G. E. Soberman, E. S. Phinney & E. P. J. van den Heuvel 

Winds Ring 
q — > q^oo q — > g^oo 



d In M T 
d In m 2 



1 9(1- a 1 9(1-5) 

a + /3 (5 

Aa-l q(l-a) 7( 5-l g( 7 5+(l-<5)) 

Ifef 3(Aa-l) 3( 7 5-l) 3 2^> 

lift 2(,4a -1) 2t ^_ 2(7* -1) 22^^ 



g m 2 M T,0 m 2 M T,0 

M T 1 22. 1 22. 

M T0 9 g 

, , , /li Aa-l r , (-«S + (l-,5) 

(L//j)o Vgo/ Vgo/ V go / V go / 



P 
Pa 



( (LM V 
V (Vm)o ,/ 

( jmulY 



Table 4. Formulae describing orbital evolution in the limits of extreme mass ratios, q = m 2 /m 1 <C 1 and q^$> 1. 
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dlnP 
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Where the relevant exponents are functions of the param- 
eters 7 and S: 



It is also instructive to examine the model in the degen- 
erate cases of e = and e = 1, where the functional forms 
change. When there is no accretion (e = 0), the standard 
C5 becomes singular, while the term 1 + eq approaches 1. 
Defining the singular part of C5: 



'sing 



lim eCs 

(3 + jS, 



the equations governing binary evolution can be rewritten, 
for the case when no material is accreted: 
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a 

P 



9 \ /I 



9o 



B 5 



9o J \ 1 + 9 

exp [(9 - qo)C s ing] , 
2(A 5 -1 



1 



9(i 



1 + 9 

exp [2C sing (q -90)] , 

X 3(^5-1) /, , 
9 \ / 1 + 9o 



9o / V 1 + 9 

exp [3C si „ s (9 - 90)] , 



1-2B5 



i-38 5 





= j5 


(B15) 


d In a 




= 


(B16) 


d In g 




7<5 2 
1 -(5 


(B17) 


<91nP 


e r 


= 1 -(5 . 


(B18) 


<91ng 



= 2(A-1) + (1-2S 5 ) T 

+ (3 + 2C si „ 9 )9 , 
= 3(A - 1) + (1 - 3B 5 ) 



+ 9 



1 + 9 



+ (5 + 3C sing )q . 



(B20) 

(B21) 

(B22) 
(B23) 
(B24) 



(B19) 



It might be worth noting that the above considerations 
are irrelevant for the pure a = 1 models, where (3, 8, and 
e vanish. There is no profound reason for this. 

In the case where all material is accreted (e = 1), 
there are seeming singularities in the cocficicnts B§ and 
C5. Proper solution of the equations of evolution in this 
case, or setting eq — > q before taking limiting values of the 
coeficients B5 and C5, shows that there is no problem at 
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all. In this case, the equations reduce to Eqs. (32), (33), 
and: 



M T 
L 

(9 In a 
d\nq 
d\nP 
dlnq 
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dm 2 
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combined 

a, 13,8 
1-a-P-S 



A Eq. (25) Eq. (B15) Eq. (B8) 
B Eq. (26) Eq. (B16) Eq. (B9) 
C Eq. (27) Eq. (B17) Eq. (BIO) 
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(Eq. (22)) 
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Table 5. This reference table is divided into three parts. First 
are the model parameter definitions. An index of equations for 
the coeficients A, B, and C, relevant to each particular model, 
follows. Last are the various formulae of orbital evolution de- 
rived in this paper. 



Fig. 8. This figure shows conservative evolution in three differ- 
ent red giant neutron star binaries, initial with m x = 1.4M , 
m g = I.OMq, and Pi = 1, 10, 100 days. Red giant coremass in- 
creases monotonically as the system evolves and has been cho- 
sen as the independent variable. The lower right hand panel 
shows the evolution of the two stars' masses, with increasing 
m a (t); the neutron star's mass (mi) increases, with comple- 
mentary a decrease in the mass of the giant (m 2 ). Orbital pe- 
riod evolution is shown in the lower left hand panel. The mass 
transfer decreases the mass of the lower mass star, and conse- 
quently widens the orbit (P increases, at constant M T and L; 
Eq. (29)). The upper right hand panel shows the evolution of 
semimajor axis (upper, solid curves) and red giant Roche ra- 
dius (lower, dashed curves). The low and high coremass ends of 
each segment in this panel are indicated by a circle and cross- 
bar, for clarity. While transferring mass, the red giant fills its 
Roche lobe, so the dashed segments shown here are consistent 
with the giant's coremass radius relation (Eq. (65)). Finally, 
the upper left hand panel shows the red giant mass loss rate. 
The dotted line labeled 'Edd' is the Eddington limit accretion 
rate of the neutron star (rh x = 4-kRxc/ ' Ktk, with an assumed 
neutron star radius Rx = 10 km and Xh = 0.70. General 
relativistic effects and the variation of neutron star mass and 
radius with accretion have been ignored.). At high coremasses, 
where the red giant's evolution is rapid, the mass loss rate can 
be far in excess of the neutron star's Eddington rate, implying 
that mass transfer is not always conservative. 



This article was processed by the author using Springer- Verlag 
ET E X A&A style file L-AA version 3. 
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Fig. 9. Mass transfer by isotropic re-emission (a = 0, j3 = 1); 
details as in Fig. fel 
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Fig. 11. Mass transfer with isotropic re-emission only so strong 
as to insure m x < m x ,Edd- For details, see Fig. pi 



Fig. 10. Mass transfer by a combination of accretion and wind 
(a = 0.8, f3 = 0). The pure wind gives similar results for initial 
periods of 10 and 100 days, but leads to instability for short ini- 
tial periods. This instability in low core mass systems is easily 
understood, since shorter period systems can evolve towards 
lower mass ratio systems, where Ci( Q = 1, /? = 0) > f eg . See 
also Fig. ^ for details. 



